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Two of the most popular approximations for the distribution of 
the steady-state waiting time, Woo, of the M/G/l queue are the so- 
called heavy-traffic approximation and heavy-tailed asymptotic, re- 
spectively. If the traffic intensity, p, is close to 1 and the processing 
times have finite variance, the heavy-traffic approximation states that 
the distribution of Woo is roughly exponential at scale 0((1 — 
while the heavy tailed asymptotic describes power law decay in the 
tail of the distribution of Woo for a fixed traffic intensity. In this pa- 
per, we assume a regularly varying processing time distribution and 
obtain a sharp threshold in terms of the tail value, or equivalently in 
terms of (1 — p), that describes the point at which the tail behavior 
transitions from the heavy-traffic regime to the heavy-tailed asymp- 
totic. We also provide new approximations that are either uniform in 
the traffic intensity, or uniform on the positive axis, that avoid the 
need to use different expressions on the two regions defined by the 
threshold. 

1. Introduction. A substantial literature has been developed over the 
last forty years that recognizes the simplifications that arise in the analysis 
of queueing systems in the presence of "heavy traffic. " The earliest such 
"heavy traffic" approximation was that obtained by Kingman (1961, 1962) 
for the steady-state waiting time for the G/G/l queue. In particular, let 
W n be the waiting time (exclusive of service) of the nth customer for a first-in 
first-out (FIFO) single-server queue (with an infinite capacity waiting room) 
fed by a renewal arrival process [with i.i.d. inter-arrival times (xn'- n > 1)] 
and an independent stream of i.i.d. processing times (V n :n > 0). If p = 
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EVi/Exi < 1, then W n => Wx> as n — > oo, where Woo can be approximated 
via 

n i\ w £Varxi + VarFi 

(L1) 2(^xi-Wi) P( } 

when /? is close to 1. Here, Exp(l) is an exponential r.v. with mean one and ~ 
denotes "has approximately the same distribution as." A precise statement 
of the limit theorem supporting the heavy traffic approximation (1.1) is given 
by (2.4) below. The term "heavy traffic" arises as a consequence of the fact 
that long queues tend to form when EV\ and E\\ are roughly balanced. 

The modern approach to justifying (1.1) involves first showing that (W n : 
n > 0) can be approximated in heavy traffic by a one-dimensional reflecting 
Brownian motion (RBM) [see, e.g., Iglehart and Whitt (1970a, 1970b)] and 
then verifying that the steady-state r.v. Woo can be approximated by that 
of the RBM [Szczotka (1990, 1999)]. Similar methods apply, in significant 
generality, to multi-station queueing networks. For example, Reiman (1984) 
proves a functional limit theorem that justifies approximating a single class 
multi-station queueing network by multi-dimensional RBM. Recent work of 
Gamarnik and Zeevi (2006) establishes the associated steady-state conver- 
gence. Harrison and Williams (1987) analyze the multi-dimensional RBM 
and show that it has exponential tails. 

On the other hand, if the processing times are heavy-tailed (e.g., regularly 
varying), there is a significant literature that establishes, for various models, 
that the associated queueing system possesses a heavy-tailed steady-state. 
A representative result of this type states that when p < 1 for the G/G/ 1 
FIFO queue described above (with regularly varying processing times), we 
have 

A f°° 

(1.2) P(Woo > x) ~ / P(Vi > y) dy 

as x — > oo, where A = 1/Exi [see, e.g., Embrechts and Veraverbeke (1982)]. 
Corresponding heavy-tailed steady-state asymptotics also exist in the con- 
text of queueing networks [see, e.g., Baccelli, Schlegel and Schmidt (1999) 
and Baccelli and Foss (2004)]. 

At first, it may seem contradictory that the heavy-traffic theory typically 
predicts exponential tails for the steady-state distribution, whereas regularly 
varying heavy-tailed asymptotics predict power-law decay in the steady- 
state tail. Of course, the key is to note that the two families of results 
involve different types of limits, one as p — > 1 (heavy traffic) and the other 
as x — > oo (heavy tails). The interesting mathematical issue here is therefore 
to send p to 1 and x — > oo simultaneously, and to determine the x- value 
(as a function of p) at which the steady-state distribution begins to "feel" 
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the presence of the heavy tails in the processing times. In particular, this 
paper develops a very explicit description, in the setting of the M/G/l queue 
(in which the arrival process is assumed Poisson), of where the transition 
from the exponential heavy-traffic approximation (1.1) to the heavy-tailed 
approximation (1.2) occurs. As a corollary to our main results (Corollary 
2.3) we find that when the processing times are regularly varying, then the 
tail probability PiWoo > x) sharply transitions at 

from the approximation (1.1) to the approximation (1.2) (where a is the 
tail index of the regularly varying V\). Roughly speaking, to the left of x* , 
(1.1) is valid whereas to the right of x* , (1.2) is appropriate. A companion 
paper [Olvera-Cravioto and Glynn (2010)] provides uniform approximations 
for P(Woo > x) in the general subexponential case, and shows how in the 
setting of Weibullian tails one can identify an intermediate zone in which 
neither the heavy-traffic asymptotic nor the heavy-tailed asymptotic hold. 

This result ties together two significant queueing theory literatures, namely 
heavy traffic theory and heavy-tailed approximations. As the first such result 
describing the transition from the heavy traffic regime to the heavy-tailed 
asymptotic, it suggests the possibility of similar such results for more com- 
plex systems and networks. Furthermore, one of our main results, Theorem 
2.1, provides an approximation for the tail probability P{Woo > x ) that is 
uniform across all values of p, and that in numerical experiments seems to 
perform very well. This new uniform approximation, which takes advantage 
of the Pollaczek-Khintchine formula for the M/G/l queue, provides a signif- 
icant numerical improvement over the existing heavy-traffic and heavy-tail 
approximations that are commonly used to approximate the tail of the r.v. 

Woo. 

2. The main results. Let (W n (p) :n > 0) be the waiting time sequence 
for an M/G/l FIFO queue that is fed by a Poisson arrival process having 
arrival rate A = p/EV\ and independent i.i.d. processing times (V n :n > 0). 
We assume throughout the remainder of this paper (unless otherwise noted) 
that V\ has a regularly varying distribution with tail index a > 2, so that 

P(Vi > x) ~ x~ a L(x) 

as x — > oo, where L(-) is slowly varying [see page 412 of Asmussen (2003)]. 

If p < 1, W n {p) => Woo(p) as n — > oo, where the Pollaczek-Khintchine 
formula [see, e.g., page 237 of Asmussen (2003)] guarantees that 

oo 

(2.1) PiWM >■) = £(1 - p)p n P{Sn > •)• 

n=0 
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Here, S n = X\ + • • • + X n (with Sq = 0), where the Xj's are i.i.d. with com- 
mon density g(-) = P{V\ > -)/EV\. The heavy-tail result (1.2) translates, in 
the M/G/l setting, into the asymptotic 

(2.2) p(w oo {p)>x)~-^P{X l >x), 

1-p 

as i-> oo. It is straightforward [see, e.g., page 404 of Asmussen (2003)] to 
show that (2.2) in turn implies that 

(2.3) P{W OD {p)>x)~-^-?—-L(x) 

1 — p a — 1 

as x — > oo . 

Turning next to the heavy traffic limit theorem for Woo{p) [due to King- 
man (1961)], its precise statement (in our M/G/l setting) is that 

(2-4) (l-p)Woo(p)=>^-Ebq)(l) 

as p /* 1, providing theoretical support for the approximation 

(2.5) P{Woo(p) > z) ~ exp(-2(l - p)EV 1 x/EV 1 2 ) 

when p is close to 1. To get a sense of the point x* = x*(p) at which the 
heavy traffic approximation (2.5) transitions into the heavy-tail approxima- 
tion (2.3), note that the point x* at which the exponential (2.5) crosses the 
power law tail (2.3) must satisfy 

EV\ 

(2.6) 2a;* (1 - p)— ± « log(l - p) + (a - 1) log**. 

-C/ V-i 

This implies that x* « «(1 -p)" 1 log((l -p)" 1 ), where re = (a- 2)EV^/(2EV 1 ). 

To make the above heuristic rigorous we look more closely at the Pollaczek- 
Khintchine formula. First we note that the heavy-tail asymptotic (2.2) can 
be obtained by simply substituting P(S n > •) by nP{X\ > •), that is, by us- 
ing the so-called subexponential asymptotic for P(S n > x). Such asymptotics 
are typically stated for fixed values of n, but can be shown to hold for n — > oo 
provided n grows slowly enough compared to x [see, e.g., Borovkov (2000); 
Rozovskh (1989)]. In other words, we can obtain the heavy-tail asymptotic 
from the first terms of (2.1), 

N(x) N(x) 

V (1 - p)p n P(S n > x) « V (1 - P )p n nP(X 1 > x) ~ -^-P(X 1 > x) 
* — ' 1 — p 

n=l n=l H 

for some appropriately defined N(x). This raises the question of whether 
we can also obtain the heavy-traffic asymptotic directly from (2.1), and the 
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answer is yes. For large n, say n > x/EXi, P(S n > x) = O(l), so by simply 
replacing P(S n > x) by one we obtain 

oo oo 

]T (i-p)p»p(s n >*)« (i-p)p n = p [x/EXl] - 

n=[x/EXi] n=[x/EXi] 

Since as p /■ 1, pl x / EX ^ ~ e -{i-p)*/EXi = e -2(i-p)EVix/EV^ we can recover 
the heavy-traffic asymptotic from the last terms of (2.1). 

This reasoning leads us to the observation that the transition of P(W OQ > 
x) occurs at the level of the partial sums P(S n > x). For the regularly varying 
case, the transition from the subexponential asymptotic nP{X\ > x — nEX\) 
to the CLT approximation 1 — — nEX\) / ' y/Vai(Xi)) [or its stable law 
counterpart when Var(Xi) = oo] occurs smoothly, which allows us to approx- 
imate the Pollaczek-Khintchine formula directly and obtain an expression 
that does not require p to be close to one. Theorem 2.1 below describes this 
(uniform in p) approximation, and Theorem 2.2 gives an equivalent formu- 
lation in terms of more familiar asymptotic expressions. As corollaries, we 
obtain the result regarding the transition from heavy-traffic to heavy-tail of 
P^Woo > x), both in terms of p as a function of x and x as a function of p. 

We also point out that similar versions of our results should also hold for 
the GI/GI/1 case. The added difficulty lies in the fact that although W OQ (p) 
still has a representation of the form 

oo 

Woo(p) =5^(1 - 0)9 n P(Y 1 + ...+Y n >x), 

n=l 

where the Y^'s i.i.d regularly varying random variables [see Asmussen (2003), 
Chapter X.9], the distribution of the Y^s and the geometric parameter 6 are 
not explicitly known. In particular, both of them depend on p, so a uniform in 
p version of Theorem 3.1 and an asymptotic expression for 6{p) are required. 
Such uniform in p results have been recently developed in Blanchet, Glynn 
and Lam (2010). Proof techniques very similar to those given here can then 
be used to obtain the GI/GI/1 equivalents of our results. 

Theorem 2.1. Suppose P(V\ > x) ~ L(x)x~ a with a > 2 and let p = 
EX 1 = EV?/{2EVi). Define p = {2 f\ (a - l))" 1 , M(x) = [{x-x p )/p\, and 

M(x) 

S(p,x)= ^{l - p)p n nP{X 1 > x - {n - l)p). 

n=l 

Then, 

P(W OQ (p)>x) 1 ^ Q 
S(p, x) + p x /v 



sup 

0<p<l 
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as x —> oo. Alternatively, 



sup 

x>0 



P(W oc (p)>x) 



S(p,x) + p x /v 



as p /* 1 . 



Theorem 2.2. Suppose P(V\ > x) ~ L(x)x a urai/i a > 2 and Zei /i 
EXi = EV^/(2EV 1 ) and j(x,p) = 1 - p x ^ - p x ^(l - p)x/p. Then, 



sup 

0<p<l 



P(W 00 (p)>x) 



{p/(l-p)h(x,p)P(X 1 >x) + p*/>* 
as x — > oo. Alternatively, 

P{W 00 {p)>x) 



sup 

a;>0 



(p/(l - p)) 7 (x, > X) + 



as p f~ 1 . 



From Theorem 2.2 we can derive the following corollary stating the dif- 
ferent regions where either the heavy-traffic approximation or the heavy-tail 
asymptotic govern the tail behavior of the steady-state waiting time. Corol- 
lary 2.3 describes the shape of the distribution of W OQ (p) for a fixed value of 
p. On the other hand, Corollary 2.4 can be of practical use in understanding 
the sensitivity of a system to the traffic intensity, since for a fixed value of 
x it tells us how P{W 00 (p) > x) changes as p gets closer to one. 

Corollary 2.3. Suppose P(V\ > x) ~ L(x)x~ a with a>2 and let k = 
(a — 2)EVi / {2EV\). Suppose that y = y(p) satisfies 

y(p) = c K (l- P r 1 log((l-py 1 ) 



for p < 1 . 

(a) 7/0 < c < I, then 



(2.7) 



sup 

0<x<y 



P(W 00 (p)>x) 



1 







exp(-2(l - p)EVix/EV?) 

as p / 1. Relation (2.7) continues to hold when c=l, provided that 
L(x)/(logx) a ~ 1 — » as x — >■ oo. 
(b) Ifc>l, then 

(2.8) "™ P{W 00 {p)>x) 



sup 

x>y 



{p/{l-p))P{X 1 >x) 



1 







as p/*l. Relation (2.8) continues to hold when c = l, provided that 
L(x)/(logx) a ~ 1 — > oo as x—>oo. 
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Corollary 2.4. Suppose P(V\ > x) ~ L(x)x a with a>2 and let k = 
(a — 2)EVi /{2EVi). Suppose that p = p{x) satisfies 

p(x) = 1 — cK,(logx)/x. 



(a) IfO < c < I, then 



(2.9) 



sup 

p<p<i 



P(W oc (p)>x) 



exp(-2(l - p)EVix/EV?) 



1 







as x — > oo. Relation (2.9) continues to hold when c= 1, provided that 
L{x) / \ogx —7-0 as x —7- oo. 
(b) If c > 1, then 



(2.10) 



sup 

0<p<p 



P(W OQ (p)>x) 



{p/{l-p))P{X 1 >x) 



1 







as x— >oo. Relation (2.10) continues to hold when c = l, provided that 
L(x)/logx — > oo as x—>oo. 



Note that Theorems 2.1 and 2.2 suggest different approximations for 
P(W OQ (p) > x). We tested both approximations and found that 

H(p,x) = S(p,x)+p x ^ 

is better than its asymptotic counterpart and performs very well for most 
values of x and p. In Section 4 we analyze how this approximation compares 
to using the simpler heavy-traffic and heavy-tail asymptotics in the regions 
where they are valid, and we give a couple of numerical examples. 

It is instructive to contrast the behavior obtained in the above regularly 
varying setting with what occurs in the light-tailed setting. Suppose, in 
particular, that E exp(9Vi) < oo for some > 0, and define 9*{p) as the root 
of pEexp(9*(p)Vi) = l. 

Theorem 2.5. Suppose that Eexp(9Vi) < oo for some 9 > 0. 

(a) Ify = y(p) = o{(l-p)~ 2 ), then 

P{W oc {p)>x) 
exp(-2(l - p)EVix/EV?) 

as p f~ 1 . 



sup 

0<x<y 
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(b) Forx>0, 
P(W 00 (p)>x(l-p)- 2 )~exp\ 

as p f~ 1 . 

(c) As P/ *l, 



A EVt EVl 



sup 

x>0 



+ X 



xEVf 
3EV? ~ 4 EV\ 



piy/^p) > x) 



exp(—9*(p)x) 



0. 



Note that in contrast to the heavy-tailed setting, the heavy traffic ap- 
proximation is now valid over a larger range, namely up to tail values of 
order o((l — p)~ 2 ). At tail values of order (1 — p)~ 2 , the third moment of V\ 
enters the asymptotic for PiW^p) > x) [see also Abate, Choudhury and 
Whitt (1995) and Blanchet and Glynn (2007)]. Finally, part (c) shows that 
the Cramer-Lundberg tail asymptotic [see, e.g., pages 365-369 of Asmussen 
(2003)] is globally valid in heavy traffic, showing the clear superiority of the 
Cramer-Lundberg asymptotic over the heavy traffic approximation when 
p is close to 1. On the other hand, for regularly varying tails, any global 
approximation to P(W OQ (p) > •) must utilize both the heavy traffic approx- 
imation and the appropriate tail asymptotic. 

We close this section with a brief discussion of how the theory described 
in this paper extends to the more general setting of geometric random sums. 
Specifically, consider the random variable 

N(p) 
i=l 

where (Yi : i > 1) is a sequence of nonnegative nonlattice i.i.d. random vari- 
ables independent of the geometric r.v. N(p) having mass function 

P(N(p) = k) = (l-p)p k ~ 1 

for k > 1 [see Kalashnikov (1997) for various applied settings in which such 
geometric random sums arise]. We assume that Y\ is regularly varying with 
finite variance, so that there exists (3 > 2 and a slowly varying function L(-) 
for which 

P(Yi > x) ~ x~ p L{x) 
as x -> oo. Put T = (P- 1)EY 1 . 



Theorem 2.6. Letp = EY\ andj(x,p) = l — (l—p) x ^ — (l—p) x ^px/p. 
Then, 

P(Z(p) > x) 



sup 

0<p<l 



(1 - p)j(x,p)P{Yi > x)/p + (1 - p) x ^ 
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as x — > oo. Alternatively, 



sup 

a;>0 



P(Z(p) > x) 



(1 -p)-f(x,p)P{Yi > x)/p + (l-p) x /» 



as p\. 0. 



Corollary 2.7. Suppose that y = y{p) satisfies 

y(p) = crp' 1 log(l/p) 



(a) IfO < c < I, then 
(2.11) 



sup 

0<x<y 



P(Z(p) > x) 



exp(— px/EYi) 

as p^O. Relation (2.11) continues to hold when c=l, provided that 
L(x)/(logx)P — > as x — >■ oo. 
(b) I/O 1, then 

(2.12) ^nz{ P )> X ) 



sup 

x>y 



P(Y 1 >x)/p 



as plO. Relation (2.12) continues to hold when c = 1, provided that 
L(x) I (log x)P — > oo as x — > oo . 

3. Proofs. In this section, we prove Theorems 2.1, 2.2 and Corollary 2.3; 
the proofs of Theorem 2.6 and Corollary 2.7 are essentially identical to those 
of Theorem 2.2 and Corollary 2.3. The proof of Corollary 2.4 is very similar 
in spirit to that of Corollary 2.3, the difference being that it follows from 
the uniform in < p < 1 statement of Theorem 2.2 instead of the uniform in 
x > 0. Theorem 2.5 follows directly from Theorem 2 in Blanchet and Glynn 
(2007). 

We now turn our attention to the proof of Theorem 2.1. Recall that 



(3.1) 



PiW^P) >x) = ]T(1 - p)p n P(S n > x) 



n=0 



where the AYs are i.i.d. with common density g(-) = P{V\ > -)/EV\ and 
S n = X\ + • • • + X n . 

Our analysis is based on the principle that we can approximate P(S n > x) 
by the heavy tail asymptotic nP{X\ > x — (n — 1)E[X]) uniformly in n 
throughout the region of large deviations of S n . Early results of this kind 
are due to Nagaev (1981), Rozovskh (1989), Mikosch and Nagaev (1998), 
Borovkov (2000), and more recently, Denisov, Dieker and Shneer (2008). The 
statement we present below is taken from Borovkov and Borovkov (2008), 
Theorems 3.4.1 and 4.4.1. 
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Theorem 3.1 (Borovkov). Let Y%,Y2,... be i.i.d. random variables hav- 
ing EY = 0, F~(t) = P(Yi > t) and F~(t) = i _/3 L(i) where L(-) is slowly vary- 
ing. Set S n = Y\ + • • • + Y n , n > 1 . 



(a) If (3 >2 and EY 2 < oo, define a(n) = ^J(f3 - 2)n\ogn. 

(b) // P G (1, 2) and F(—t) < cF{t) for t > and some constant c> 0, define 
a(n) =F~\l/n). 

Then, there exists a function p(t) J, as t f oo such that 

P(S n > y) 



sup 

y>ta(n) 



nP(Y 1 > y) 



1 



<¥>(*) 



uniformly in n. 



Below we give an application of Borovkov's result to our particular setting. 



Lemma 3.2. Let X\,Xi, ■ ■ ■ be i.i.d. nonnegative random variables with 
a = E[X] < oo, and P(X\ > t) = t~ a+1 L(t) where L(-) is slowly varying and 

a>2. Set S n = X 1 ^ \-X n , n> 1. For any (2 A (a - l))" 1 <7 < 1 define 

M 7 (x) = |_(x — x 7 )//ij . Then, there exists a function ip(t) ^ as t \ oo such 
that 

P(S n > X) 



sup 

l<n<M 7 (z) 



nP(X 1 > x - (n - l)/x) 



1 



< (f(x) 



Proof. Suppose first that a > 3 and let cr(n) = y/(a — 2)nlogn. Since 
P{S n >x) _ P(S*>x-nu) 

nP{X\ > x — (n — l)u) nP(Y\ > x — nu) ' 

where Yi = Xi — u and = Y± + ■ — \- Y n . Then the result will follow from 
Theorem 3.1(a) once we show that (x — n^)/a{n) — > oo uniformly for 1 < 
n < M 7 (x). To see this simply note that 

x — nu x — M 7 (x)/i / Jl x 7_1 / 2 
a(n) ~ a(M 7 (x)) V « ~ 2 ' v 7 * !^ ' 

Since 7 > 1/2, the above converges to infinity. 

Suppose now that a G (2, 3) and note that P(Y\ < —t) = for t > fi. 
Note also that since F(t) = P{Yi > t) is regularly varying with index a — 1, 
then <r(n) = F (1/n) = n 1 ^ " 1 ^ L(n) for some slowly varying function Z(-) 
[see Bingham, Goldie and Teugels (1987)]. Then the result will follow from 
Theorem 3.1(b) once we show that (x — nfj,)/o~(n) — > 00 uniformly for 1 < 
n < M~(x). To see this note that 

x — n/j, x — Mj(x)(i x 7 j; 7-1 /^ -1 ) 
a{n) ~ <r(M 7 (x)) ~ a(x/fi) ~ /x-V(°-i)Z( x ) ' 
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and since 7 > l/(a — 1) the above converges to infinity. 

The case a = 3 is rather technical and does not provide additional insights. 
We refer the reader to the internet supplement Olvera-Cravioto, Blanchet 
and Glynn (2010) for the details. □ 



We now give a lemma that will allow us to transform the statements of 
the main results from being uniform in < p < 1 to being uniform in x > 0, 
under the limiting regimes x — > 00 and p 1 , respectively. 



Lemma 3.3. Suppose that 

P(W 00 (p)>x) 



sup 

0<p<l 



A(p,x) 



as x—^oo, where A(p,x) satisfies 



sup \A(p, x) — lj — > 

0<x<{l-p)-n 



as p /* 1 for some < rj < 1 . Then, 

P{W 00 {p)>x) 



sup 

x>0 



A(p,x) 



as p /~ 1 . 



Proof. We argue by contradiction. Suppose that there exists an e > 
and a function x : (0, 1) — > (0, 00) such that x{<f>) > (1 — 4>)~ v and 



A{cj>,x{<t>)) 



> e 



for all < 4>< 1. Then, 

P{W 00 {p)>x{<j>)) 



sup 

0<p<l 



1 



> 



P(Woo($>s(0)) 



A(p,x(4>)) ~ A(4>,x(<t>)) 

But this cannot be since by assumption, 

PiWM > x(4>)) 



> e. 



It follows that 



lim sup 

</V*l0<p<l 



sup 



A{p,x{4>)) 
P{W 00 {p)>x) 



0. 



A(p,x) 
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as p /* 1. For < x < (1 — p) v note that 

lim sup — r 1 

P/" l 0<x<{l- P )-n A(p,x) 

, v \P(W OQ (p)>x)-l\ 

< lim sup j7 r h lim sup — ; r — 

P/ il 0<x<(l-p)-v A(p,x) P^Q<x<{l- P )-n A(p,x) 

= lim sup \P(Woo(p) >x)-l\. 

P/* 1 0<x<(l- p)-n 

The last limit is zero by the standard heavy traffic limit. □ 

Throughout the rest of this section let p = EX\, (3 = (2 A (a — 1)) , 

M(x) = [(x — x@)/p\ , and 

M{x) 

S(p, x)=^2(l- p)p n nP(X l >x-(n- l)p). 

n=l 

Proof of Theorem 2.1. We will prove the uniform in < p < 1 
asymptotic, since the statement regarding the uniformity in x > will follow 
from Lemma 3.3 by noting that 

(l-p)M{x){M(x) + l) 
sup b (p, x) < sup 

0<x<(l-p)" 1 /4 0<x<(l-p)- 1 /4 2 



< 



P 2 



1/2 



which clearly converges to zero. Throughout the proof C > is a generic 
constant. 

Fix < 7 < 1 A P(a - 1) and define M y (x) = [(a: - x 7 )//uj. Then, by 
Lemma 3.2, there exists a function — > as i — > oo such that 



sup 

l<n<M 7 0) 



P(5 n > X) 



nP(Xi > x - (n - 



1 



< ^i(x). 



By (3.1) we have 

\P(W oc (p)>x)-S(p,x)-p x ^\ 



M 7 (x) 

< ipxix) £ (1 - p)p n nP(Xi > x - (n - l)/z) 

n=l 
M(x) 

+ £ (1 - p)p n nP(X 1 >x-(n- l)p) 

n=M 1 (x)+l 
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+ ]T (i-p) P n p(s n > 

n=M 1 (x)+l 



+ 



(1 - p)p n P{Sn >X)~ 

n=lx/fi\+l 



Clearly, 



M~ l (x) 



(3.2) tpxix) p)p n nP{X l >x-(n- l)p) < ipx (x)S(p, 



Fix < e < min{a — 2, (fi(a — 1) — 7)//?}. The second term is bounded by 

M(x) 

{I- p)p n nP{X 1 >x-{n-l)p) 

n=M 1 {x)+l 

M(x) 

<Cx(l-p) P n {x-{n-l)p)- a+1+e . 

n=Mj(x)+l 



Since g(n) = p n (x — (n — l)p) 



-a+l+e ; 



is convex m n, 



x 



M{x) 

(1-p) 9{n) 

n=M j(x)+l 

< x(l - p){M(x) - M 7 (x))max{si(M 7 (:r) + l),g{M(x))} 

(3.3) 

< C(l - p)x 1+7 nlax {p(*-*' r )/M a; -7(a-l-e) } p(*-* /, )/M-l :c -/'(a-l-e) J. 

< C(l - p ) a .l+7-^(«-l-e) p (*-xT)/M j 

where our choice of e guarantees that 7 — /3(a — 1 — e) < 0. Also, we have 

L^/mJ 

Y (1 - p)p n P{S n >X)< p ^7^)+l(l _ p L*/Mj--Mnf(*)) 
n=M 7 (i)+l 

(3.4) 

<C7,( 3S -* r )/'V|logp|. 
To derive the last bound let K(x) = [(a; + x 7 )//ij . Then, 



£ (l-p)p n P(S n >x)-p^ 

n=[x/fi\+l 
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= ^//i_ p L*/Mj+l+ {l-p)p n P{S n <x) 
n=[x//fj+l 
K{x) oo 
<f?/»{l-p)+ ^~P)p n + £ (l-p)p n P(Sn<x). 

n=[x/fi\+l n=K(x)+l 

It is easy to check that 

K(x) 

£ (l-p)p n <Cp*/W\logp\. 

n=\x/ 

For the tail of the sum let Yi = p — Xi and S*=Y\ + -- - + Y n . Let b n be 
the scaling for which Z n = S*/b n =>■ Z, where Z is a stable random variable. 
Note that b n = n^Lo(n) for some slowly varying Lq(-). It follows that for all 
n > K(x), 



P(S n <x) = P[Z n > 



Tip — X 



<P[Z n > 



where 



(K(x) + l)p-x 



> 



prx 1 



(K(x) + l)p- x 

r7~P 



> 



cx ' 



{x + x~/)PL {(x-x7)/p) ~ L (x) 



for some constant c > 0. It follows that 

sup P(S n <x)< sup P(Z n > cx 1 ~^/Lq(x)) < if2(x) 

n>K{x) n>K(x) 

for some <p2(t) — > as t — > oo. Hence, 

oo oo 

£ (i-p)p n P(s n <x)<Mx) £ (i- P )p n <Mx)p x/l1 . 

n=K{x)+l n=K(x)+l 



We thus have that 



(3.5) 



£ (l-p)p n P(S n >x)-p^ 

n= [x/fj.\+l 



<^ 2 {x) P x ^ + C P x ^\\ogp\x\ 
Combining (3.2)-(3.5) gives 

\P(W oc (p)>x)-S(p,x)-p x ^\ 

< Mx)S( P ,x) + Mx)p x/ ^ + Cp^-^\logp\x v , 
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where v = max{l + 7 — f3(a — 1 — e), 7} € (0, 1). It only remains to show that 
p(x-x~i) 1 >| i gp| x « = o(S(p,x) + p x /^) uniformly in < p < 1. To see this let 
p(x) = 1 — (x^logx) -1 , then 

snp ^7 7) t |l0g/ ; |X " < sup e l^^|logp|x« 

p(x)<p<l S(p,x)+p x IV p(x)<p<l 

= e n°£p(x)\*y^i ogp ( x )\ x v 
c 

log 2; 

and since 5(p, x) > P(Xi > x) J (1 ~ P)P n , 

sup 



o<p< P (x) S(p,x) + p x /» 

p(a=-^)/M|lo g p|^ 



- 0< ;5 I) i(i)x-«+v(i-/ ) W''j) 

<C sup p( a - a;7 )/^ 1 |logp|x ,;+Q - 1+£ 

0<p<p(z) 

< C sup expf— (— — 1 J t + logi + (v + a — 1 + e) log a; J 

= Cexpf — ^ C 1 - - loglogx + (a - 1 + e)logx 

\ piogx \ X 1 1 X J 

-)-0. □ 
We now prove Theorem 2.2. 

Proof of Theorem 2.2. Again, we only prove the statement regard- 
ing the uniformity in < p < 1 , since the statement for x > follows from 
Lemma 3.3 and the observation that, as p /* 1, 

P ^- p ^-^-P>p^\ P{Xl 



sup — : — — p""* \F{X 1 >x) 

0<x<(l-p)- 1 /4 1 - P V P / 

P fn ,X (l-p)x i 
< sup - J — [\\ogp\ t -!—p x ll 1 



0<X<(l-p)-V4 1 - P V M M 



Ml-P) 5 / 4 

o((i-p) 1/2 ). 



(|logp|-(l-p)e" |los ' ,|(1 " p) 1/4/ ^) 
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By Theorem 2.1 we only need to show that 

S(p,x)-(p/(l-p)h(x,p)P(X 1 >x) 



sup 

0<p<l 



(p/(l - p)) 7 (x, p)P(X 1 >x) + pFl» 
as x — > oo. We start by noting that 



S(p,x) 



-j(x,p)P(X 1 >x) 



1-p 

M{x) 

(1 - p)p n n{P{X l >x-(n- l)p) - P(X 1 > x)) 



n=l 



/M(x) \ 

+ P{X 1 >x)[ Yj (! " P)P Hn ~ TTp- l{x > p) ) 
Then, since j^j{x,p) > En=f' J ~ 1 ( 1 ~ P)p n ^ 



l-p 



j(x,p)P(X 1 >x) 



M(x) 



< P(X, >*)£(!- pKn( P ( X '>^-("r'W - 

n=l ^ 



> x) 



+ p(Xi > x) \jz-pi(^ p) - E o- - p)p n A 

The second term can be bounded as follows: 

/ M(x) \ 

P(Xl > x) Y=-l{x, p) - (1 - 

V n=l / 



(3.6) 



< -^-P(Zi > x)(p M ^ - p^) ( 1 + 



< 



P_ 
l-p 

C 
l-p 



(l-p)x 



/' 



P{X l > x){p^ x - x ^^ - p x ^ +1 )(l + (1 - p)x), 



where C > is a generic constant. Fix < e < j3{a — 2)/(a — 2 + /3) and set 
JV(x) = L(l - e)x/p\ ■ Then, for 1 < n < N(x), 



P(Xi >x-(n- l)p) ( x - (n-l)p 



P(Xx > x) 



x 



-a+l 



sup 

l<n<N(x) 



L(x — (n — l)p) 



< l + (a-l)e 



-a-2 



(n - l)p 



(1 + Vi(x)), 
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where <fi(x) = sup e<t<1 ffiv — 1 — > by properties of slowly varying func- 
tions. Therefore, for 1 < n < N(x), 

P{Xx >x-(n- l)/i) 



P(Xi > x) 



1<C 



n — 1 



+ <pi(x) . 



It follows that 



n=l \ \ / / 

< C7P(Xx > x) - (1 - p)p»n(n - 1) + ¥>i 0*0 ^7^, p) 

\' r n=l ^ / 

(3.7) < CP(X 1 > x) ( .- (1 - p^*) - (1 - p)N(x)p"to) 



<^-P{X 1 >x) 1 (x,p) 



1 



.ac(l-p) 

For the terms iV(cc) <n< M{x) we have 
AT (a;) 



+ <pi(x)- 70e>/0 

1-p 



p(Xl>x) S a- ^'^V^ -i) 



n=iV(a;)+l 



Af(x) 



(3.8) 



< Cx(l - p)p N ^ +1 P{Xx>x-{n-l)p) 

n=N(x)+l 



< Cx(l - p)p( 1 -c)x/» 



x-fiN(x) 



P(Xi > t) dt 

' x—/iM(x) 

< Cx(l - p)p^- £)x ^(x - pM(x))P(X 1 > x - pM(x)) 

< C(l - p)p( 1 - £)x ^x 1+ Pp(X 1 > x 13 ), 

where for the third inequality we used Proposition 1.5.10 in Bingham, Goldie 
and Teugels (1987). Combining (3.6)-(3.8) gives 

S(p,x)--^-(l-p^)P(X 1 >x) 
1-p 

< -r^—PiXi > x){p^ x - x ^^ - p x ^ +1 )(l + (1 - p)x) 
1-p 
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+ 



Cp 
1-p 



P(Xi>x)rf(x,p) 



1 



x(l-p) 



+ <Pi(x) 



+ C(l - p)p {l - £)x/ ^x 1+p P{X 1 > 



a 



Let A(p,x) = ^-j(x,p)P(Xi > x) + p x ^ and define p(x) = 1 — cp( 

2)logx/x, with ^(a^)*^-^ < c < 1. Note that j(x,p) ~ 1 as x — > oo uni- 
formly for < p < p(x). Then, 



1 



sup 



o<p<p(x) A(p,x) 



S(p,x) 



P 



1-p 



(1 - p x ^)P(X 1 > 



< C sup 

0<p<p(x) 



_J_ {p (x- x ^-i _ pX/f , ){1 + p)x) + __L_ 

J{x,p) x(l-p) 

+ + ^(x^P^x) P J 



< ci + (i _ p(x)) x ) 



+ 



x(l — p{x)) 



+ ^^ x^{l- P {x)fp{xt-^ 



<c 



logx 1 . . P(Xi > x 



(logx) 5 



P(Xi > X) s l-^+c(«-2)(l-e) 



The first three terms in the expression above clearly converge to zero. To 
see that the fourth one does as well use Potter's theorem [Bingham, Goldie 
and Teugels (1987), page 25] to obtain 

P(X 1 >x p ) (logx) 2 



P(X 1 > x 



x l-l3+c(a-2)(l-e) 

(logx) 2 



x l-P+c(ct-2){\-e) 

for some constant A e > 1. Our choice of e and c guarantees that 1 — (3 + 
c(q - 2)(1 - e) > (1 - 0)[a - 1 + e). 

To analyze the supremum over p(x) < p < 1 we first note that 



Lx/ M J+l 
n=l 



(1-P) 



Then, 



sup 



p(x)< P <i A(p, : 



S(p,x) 



1-p 



(1 - p x ^)P(X 1 > 
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<C sup {^(^-p)(l + (l- pW 

p(x)<p<\V l ~P 



P1 {x lP )P{X l >x) ( 1 

+ (1 - p)p- £x/lM x 1+l3 P(X 1 > x 13 ] 
<C sup J P ^ Xl > ^ | log p\ x 13 log x + p- £X ^x p log xP(Xi > 

p(x)<p<l [ 1 _ P 



+ > X) £ P n - x/ "n(l + x(l-p)) 



n=l 



< C{ P{X 1 



> x)^^ log x + p(x) ex lv x P log X P(X 1 > x 13 ) 
+ P{X l > x) Y p{x) n ~ x/fl (l + log x) \ . 

n=l J 

The first term clearly converges to zero. To see that the second and third 
terms converge to zero as well note that 

p{x)- £X ^xHogxP{X 1 > x 13 ) < Cx e < a - 2 ^-^ a - x - e hogx 

< Cx e(a-2)-P(a-2-E) bg x 

and 

\x/p~\+l \\x/ii\+l\ 

P(X 1>X ) p(x) n - x/ "logx<Cx- a+£ \ P[X) — Uogx 

< Cx~ a+£+l+c ( a ~ 2 ^ 

< Cx~ l+£ . 

Our choice of e guarantees that both expressions above converge to zero. 
This completes the proof. □ 

We end this section with the proof of Corollary 2.3. 



Proof of Corollary 2.3. Let 

y = y(p) = cp(a - 2)(1 - p)~ l log((l - p)" 1 ). 
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We start with the proof of part (a). We need to verify that for < c < 
(p/(l-p)) 1 (x,p)P(X 1 >x) + p x ^ 



sup 

0<x<y 



as p f - 1. Note that 



exp(-(l - p)x/p) 



sup 

0<x<y 



(p/(l-p)) 7 (x,p)P(X l >x) + p*/» 



(3.9) 



(3.10) 



< sup 

0<x<y 



exp(-(l -p)x/p) 
' x\ogp (1 — p)x 



exp 



+ 



+ sup 



P 

fTf(x,p) 



/' 



1 



0<x<y (1 - P) 

We can bound (3.9) as follows: 

' x\ogp (1 — p)x 



sup 

0<x<y 



exp 



+ 



P(Xi > x) exp((l - p)x/p). 

x\\ogp+l- p\ 



1 



< sup 

0<x<y 



P 



p p 

<Cy(l-pf 

<cr l \ogt, 

where t = (1 — p)~ l . Also, note that since p x l ^ > 1 — | log p\x/p and | log 
l- p + 0({l-p) 2 ) asp SI, 

1 (x,p) = l-p x ^-p^(l-p)x/p 

< (I logp| - (1 - p))x/p + I logp|(l - p){x/pf 

<C(l-p) 2 x 2 . 

Then (3.10) is bounded by 

pj(x,p) 



sup 

0<x<y U — PJ 



P(Xi > x) exp((l - p)a;///) 
7 (x,p) 



< C sup 

0<X<(l-p)- 1 /4 1-P 
1 

+ sup 

(l-p)- 1 /^^ J- — p 

< sup C(l - p)x 2 

0<x<(l-p)- 1 /4 



P(Xi > x) exp((l - p)x//u) 



+ sup 



exp((l — p)xj p — (a — 1) logx) 
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<C(l-p) 1/2 + 



1-p 



exp((l - p)y/fi - (a - 1) log y) 



< Ct~ 1/2 + CL(t log t) exp(- (1 - c) (a - 2) log £ - (a - 1) log log £) . 

Clearly, if < c < 1, then the two expressions above converge to zero as 
t —7- oo. If c = 1 and L(x)/(logx) Q_1 — > as x — > oo, then 

(p/(l - p)) 7 (x, > x) + fl* 



lim sup 



exp(-(l - p)x/p) 



< C lim L(t logi) exp(— (a — 1) loglogt) 

t— >oo 



C lim 



L(ilogi) 



log(tlogt) 
logt 



a-1 



t^oo (log^logi)) " 1 
We now move to part (b) . We need to verify that for c > 1 
(p/(l - p)) 7 (x, p)P(X 1 > x) + 



sup 

x>y 



(p/(l-p))P(X 1 >x) 
as p f - 1 . Note that 

(p/(l - p)) 7 (x, p)P(Xi > x) + p*/" 



sup 



1 



(p/(l-p))P(X 1 >x) 

< sup |7(x, p) - 1| + sup ——— 

x>y x>y P-T l A l > X) 

< p v/ »{l + (1 - p)y/p) + Csup -77— r exp 



L(x) 



--(1 - p) + (a - l)logx 



< Cp^(l - p)y + C^—£ exp ( - ^(1 - p) + (a - 1) logy 



<Ct" c(a - 2) log£ 
C 



+ 



exp(— (c — l)(a — 2) logt + (a — 1) loglogi), 



L{t\ogt) 

where t = (1 — p) _1 . Clearly, if c > 1 the above converges to zero. If c = 1 
and L(x)/(logx) a_1 — > 00 as x— t-oo, then 

(p/(l-p)) 7 (x,p)P(X 1 >x)+p^ 



lim sup 

P/' 1 0<x<y 



< C lim 



1 



t->oo L(t log t) 



exp(-(l -p)x/p) 

exp((a — 1) log log t) 
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4. Numerical approximations. Theorems 2.1 and 2.2 suggest approxi- 
mating P(Woo{p) > x) either with 

M{x) 

H(p, x) = S{p, x) + p x ^ = £ (1 - p)p n nP(X 1 > x - (n - + p x ^ 

n=l 

or with 

J(p, x) 4 -P-jip, x)P(X 1 >x) + pP/n, 
1-p 

respectively. 

We compared both approximations to simulated values of P(W 00 (p) > x) 
and found that H (p, x) tends to be better than J(p, x) and seems to perform 
very well across all values of x for different choices of p. This is not surprising 
given that H(p,x) more closely resembles the Pollaczek-Khintchine formula 
than J(p,x). 

When a 2 = Var(Xi) < oo, the central limit theorem can be used to ap- 
proximate the tail of the Pollaczek-Khintchine formula in a way that a 2 
is incorporated into the approximation. The term p x ^ appearing in the 
definitions of H(p,x) and J(p,x) can be replaced by 

oo 

T(p,x) ^ ]T(1 - p)p n (l - <&((* - np)/V^)), 

n=l 

which can alternatively be written as T(p,x) = E[p M ( x,z ^], where Z ~ A^O, 1) 
and M(x,z) = [(^/x/p + (az) 2 /(2p) 2 — (az) / (2p)) 2 \. We do not give proofs 
here, but it can be shown that provided a 2 < oo, Theorems 2.1 and 2.2 
continue to hold with p x ^ replaced by T(p,x). This is relevant from the 
numerical standpoint since the resulting approximations tend to perform 
better than those with the simpler p x ^ . 

We plotted approximation H (p, x) against simulated values of P(W 00 (p) > 
x). Figures 1 and 2 correspond to queues having Pareto integrated tail dis- 
tribution, that is, P(Xi > x) = x~ a+1 for x > 1. For comparison purposes 
we also plotted the heavy-traffic approximation, 

Heavy-Traffic = exp(— (1 — p)x/p), 

and the heavy-tail asymptotic, 

Heavy-Tail = —?—P(Xi > x). 
1-p 

The vertical line corresponds to the value 

x09) = /i(a - 2)(1 - p)- 1 log((l - p)- 1 ). 

The simulated values of P(Woo(p) > x) were obtained using the condi- 
tional Monte Carlo algorithm from Asmussen and Kroese (2006), and each 
point was estimated using enough simulation runs to obtain a relative error 
of at most 0.05 with approximately 99% confidence. 
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Fig. 1. Pareto integrated tail with p = 0.95 and a — 3.1. 
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